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ABSTRACT 


Soft coral Dendronephthya gigantea (Verrill, 1864) is a conspicuous species dominating shallow sea waters 
of Jejudo Island, Korea. Recently its whole mitochondrial genome sequencing was completed by us and the 
sequence information provided an opportunity to test the age of Octocorallia and time of evolutionary separa- 
tion between some representative orders of the subclass Octocorallia. Molecular phylogenetic analyses based 
on 13 mitochondrial protein encoding genes revealed a polyphyletic relationship among octocorallians repre- 
senting two orders (Alcyonacea and Gorgonacea) and four families (Alcyoniidae, Nephtheidae, Briareidae, 
and Gorgoniidae). Estimates of divergence times among octocorallians indicate that the first splitting might 
occur around end of or after Cretaceous period (50-79 million years ago (Ma)). The age is relatively young 
compared to the long history of stony sea corals (>240 Ma). Taken together our result suggests a possible 
relatively recent radiating evolution at least in the order Alcyonacea and Gorgonacea. Molecular dating and 
phylogenetic analysis based on much broader taxon sampling and many genes might give an insight into this 


interesting hypothesis. 
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INTRODUCTION 


The soft coral D. gigantea (Verrill, 1864) is a dominant spe- 
cies around Jejudo Island, Korea. This species belongs to 
the subclass Octocorallia. Understanding of evolutionary 
relationships among octocorallians remained poorly until 
near past and still unresolved in most basal nodes of phy- 
logeny. A few previous phylogenetic studies identified dis- 
tinct clades within the taxonomic group, but some clades 
did not agree with the conventional taxonomy (Won et al., 
2001). Won et al. (2001) found a lack of resolving power in 
18S rDNA data within the subclass Anthozoa. Recently 
McFadden et al. (2006) assessed evolutionary relationships 
of octocorallians using two mitochondrial genes (ND2 and 
MSH1) from the largest taxon sampling (103 genera and 28 
families) to date. This comprehensive study clearly showed 
an inconsistency between taxonomy and phylogeny, alth- 
ough basal relationships among major clades of octoco- 
rallians could not be resolved (McFadden et al., 2006). 
Largely difficulties of morphological classification in this 
group lie behind this discrepancy between taxonomy and 
phylogeny. On the other hand, application of molecular 
markers for phylogenetic analyses has been delayed in the 
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Octocorallia due to dearth of appropriate markers. In parti- 
cular evolution of mitochondrial DNA sequences in Antho- 
zoa turned out to be extremely slow (Shearer et al., 2002). 
Thus mitochondrial genes lack sufficient resolution for dis- 
criminating species within genera (Sanchez et al., 2003; 
Cairns and Bayer, 2005; Wirshing et al., 2005). Nuclear 
18S rDNAs did not provide resolution for deeper nodes with- 
in trees (Won et al., 2001). Recently Medina et al. (2006) 
applied whole mitochondrial genome to test evolutionary 
relationships among stony corals. Based on the phylogeny 
reconstructed from 13 different protein encoding genes, 
molecular dating of diverging ages at basal nodes was achi- 
eved. Because of large amounts of informative site in the 
long DNA sequences, the evolutionary relationship at the 
ordinal level between the order Scleractinia and Coralli- 
morpharia could be clearly shown. However, more recently 
other research group who used small set of genes but larger 
taxon sampling found a contrasting result (Fukami et al., 
2008). This conflict leaves further evaluation. 

Our purpose in this study was to construct a well resolved 
phylogeny among octocorallians including the dominant soft 
coral D. gigantea and to estimate divergence times among 
them. We updated the earlier molecular data set made of 
most protein encoding genes of mitochondrial genome (Me- 
dina et al., 2006) with newly sequenced mitochondrial ge- 
nome of D. gigantea. This new data set consists of main 
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Table 1. Species and their mitochondrial genes used in the study 


Species GenBank Notes* Referencest 
Octocorallia 
Dendronephthya gigantea FJ372991 = (Kim, 2008) 
Pseudopterogorgia bipinnata DQ640646 — (Medina et al., 2006) 
Sarcophyton glaucum AF064823, AF063191 No sequences for ND3 and ND4 (Pont-Kingdon et al., 1998) 
Briareum asbestinum DQ640649 - (Medina et al., 2006) 


Corallimorpharia 
AB441341 for CYTB 


Corynactis californica AB441256 for COI 


Ricordea florida DQ640648 

Discosoma sp. DQ643965 

Rhodactis sp. DQ640647 
Scleractinia “Complex” clade 

Porites porites DQ643837 

Siderastrea radians DQ643838 

Acropora tenuis NC_003522 

Pavona clavus DQ643836 

Agaricia humilis DQ643831 
Scleractinia “Robust” clade 

Astrangia sp. DQ643832 

Montastraea faveolata NC_007226 

Mussa angulosa DQ643834 

Colpophyllia natans DQ643833 
Actiniaria 

Nematostella sp. DQ643835 

Metridium senile NC_000933 
Zoanthiniaria 

Palythoa sp. DQ640650 


No sequences for CYTB, ND5, ND1, 


Only available 
for CYTB and COI 
— (Medina et al., 2006) 
— (Medina et al., 2006) 
— (Medina et al., 2006) 


(Fukami et al., 2008) 


— (Medina et al., 2006) 
— (Medina et al., 2006) 
= (van Oppen et al., 2002) 
— (Medina et al., 2006) 
— (Medina et al., 2006) 


— (Medina et al., 2006) 
= (Fukami and Knowlton, 2005) 
— (Medina et al., 2006) 
— (Medina et al., 2006) 


— (Medina et al., 2006) 
= (Beagley et al., 1995) 


ND3, COII, ND4, and ND6 (Meēdina:etal; 2006) 


*Unless noted otherwise in the column, all the sequences from the thirteen protein encoding genes were included for phylogenetic analyses: three 
cytochrome c oxidases (COI, COII, COIII), cytochrome b (CYTB), and 7 NADH dehydrogenases (ND1, 2, 3, 4, 4L, 5, 6), and ATP synthases (ATP6 and 


ATP8). 


tOriginal publications that reported and deposited the mitochondrial genome sequences. 


taxonomic representatives of anthozoan groups which are 
closely related to octocorallians. 


MATERIALS AND METHODS 


Phylogenetic Analysis 


We conducted phylogenetic analyses using thirteen mito- 
chondrial protein encoding genes from 20 taxa including D. 
gigantea. Recently we completed mitochondrial genome 
sequencing of D. gigantea and deposited the sequence infor- 
mation to GenBank (accession no. FJ372991) (Kim, 2008). 
DNA sequences of the other 19 taxa representing Hexaco- 
rallia (Scleractinia, Corallimorpharia, Actiniaria, and Zoan- 
thiniaria) and Octocorallia were collected from GenBank 
(Table 1). The single largest data set on whole mitochondri- 
al genomes of anthozoans was published by Medina et al. 
(2006). We included those sequences in the present phylo- 
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genetic analyses. Alignments of the sequences were gener- 
ated automatically and then edited manually based on amino 
acid translation in MEGA v.3.1 (Kumar et al., 2004). Ambi- 
guous sites in homology due to indels were manually exclud- 
ed and sites with gaps were included unless over half of the 
taxa are missed in the sites. A single frame of amino acid 
translation was applied to the concatenated data so that a 
data partition according to codon positions is easily made 
for Bayesian analysis. The final data set resulted in a con- 
catenated sequence of 11,946 nucleotides (Table 2). To infer 
evolutionary relationships among the taxa, maximum likeli- 
hood (ML) trees and Bayesian trees were searched by PAUP 
v. 4.10b (Swofford, 2002) and MrBayes v. 3.12 (Huelsen- 
beck and Ronquist, 2001), respectively. The best-fit model 
of nucleotide substitution for ML search was chosen based 
on Akaike information criterion (AIC) using Modeltest pro- 
gram v. 3.7 (Posada and Crandall, 1998). According to the 
model testing, a general time reversible model (GTR) with a 
proportion of invariable sites (I) and gamma (T`) distributed 
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Table 2. Thirteen protein encoding genes examined and their 
final length of DNA sequences for phylogenetic analyses 


Gene Length (n.t.) 
COI 1,542 
COII 744 
COIII 786 
CYTB 1,521 
ND1 966 
ND2 1,095 
ND3 345 
ND4 1,431 
ND4L 300 
NDS 1,761 
ND6 558 
ATP6 693 
ATP8 204 
All 11,946 


rate variation among sites (GTR+I+ T model) was chosen 
for our data set. ML search was run under the following set- 
ting: tree bisection and reconnection (TBR) branch swap- 
ping with 100 random addition and 100 bootstrap replicates. 
Bayesian analyses were conducted under general time rever- 
sible model (GTR) with variable rate partitions for each 
codon position with gamma distributed rate variation among 
sites (GTR+SS +T model) as follows: Markov chain Monte 
Carlo (MCMC) with 2 million generations, 6 parallel chains, 
sampling every 1000 iterations, and burn-in 500. We used 
the last 1,500 sampled trees to estimate Bayesian posterior 
probabilities (BPP). 


Estimation of Divergence Time 


In order to estimates divergence times of each node of the 
phylogenetic tree, we used R8S v. 1.7 program by Sander- 
son (2002). The ML tree with branch lengths was used in 
the analysis of R8S. Because there is no known appropriate 
outgroup which has sequence information of all 13 protein 
encoding genes examined in the ingroup, we chose midpoint 
rooting method. A round of cross validation procedure test- 
ed log smoothing parameter value from 0 to 10, and found 1 
the best parameter value for the data set. The value is su- 
fficiently low to allow rate heterogeneity across branches. 
To calibrate divergence time we followed Medina et al.’s 
schemes (see Fig. 1, 2006): one from “robust” clade and the 
others from “complex” clade which represent two highly 
divergent clades of the order Scleractinia (Romano and 
Palumbi, 1997a, b; Chen et al., 2002). The former is genus 
Astrangia (70 Ma: an7 node in Fig. 1) and the latter are 
genera Pavona (33 Ma: an8 node) and Acropora (55 Ma: 
an9 node) (Fig. 1). Additionally a fossil record of the first 
appearance of Scleractinia at 240 Ma was set as a lower date 
boundary for the node an16 (Fig. 1). Labeling of the nodes 
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Corallimorpharia 


Scleractinia 
“Complex/Long” clade 


Astrangia 
Montastraea 

Scleractinia 

Mussa “Robust/Short™ clade 

an5 

Colpophyflia 

Nematostella 
Actiniaria 
Metridium 


cocnnsanenanenasessan ==- Palythoa ]Zoanthiniaria 


Briareum (Gorgonacea) 


Sarcophyton 

(Alcyonacea) | Octocorallia 
| ergosacea) 

Dendronephthya 

-e {Alcyonacea) 

Fig. 1. Maximum likelihood (ML) tree of some representative 
anthozoans based on concatenated thirteen mitochondrial pro- 
tein encoding genes (11,946 n.t.). The root of the tree indicates 
a midpoint of the longest branche. Lengths of branches are pro- 
portionally drawn to their substitution rates. Numbers above or 
below branches indicate bootstrapping values of ML and Baye- 
sian posterior probabilities (BPP). If both values are equal, they 
are represented as a single value. Dashed lines indicate taxa 
that have incomplete protein encoding genes. Names of each 
node are also placed between two bifurcating branches from 
“ani” to “ani6” in series. Fixed divergence ages based on ear- 
liest fossil records are shown at nodes indicated by open circle 
for Astrangia (70 Ma) of the “robust” clade and by filled circles 
for Pavona (33 Ma) and Acropora (55 Ma) of the “complex” clade 
of Scleractinia. 


included in the R8S analysis is represented in Fig. 1 and 
Table 3. 


RESULTS AND DISCUSSION 


We conducted maximum likelihood and Bayesian phylo- 
genetic analyses involving the newly sequenced temperate 
soft coral species, D. gigantea (Alcyonacea: Nephtheidae), 
and other 19 anthozoan taxa representing Hexacorallia (Scl- 
eractinia, Corallimorpharia, Actiniaria, and Zoanthiniaria) 
and Octocorallia (Fig. 1). The phylogenetic analyses result- 
ed in the same tree topology between ML tree (Fig. 1) and 
Bayesian tree (data not shown). Almost all nodes in the ML 
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Table 3. Estimated divergence times for each node in millions 
of years 


Node Age 
ani 27-43 
an2 38-60 
an3 50-79 
an4 238-380 
an5 6-9 
an6 19-30 
an7 70* 
an8 33* 
an9 55* 
an10 70-108 
anil 94-148 
ani2 16-26 
ani3 71-114 
ani4 105-168 
ani5 192-305 
ani6 240-450 


Ages of root and the next two deeper nodes from which the order Actinia- 
ria and Zoanthiniaria descend are not shown because of inaccuracy of 
method in determining age of root and a low ML bootstrapping support 
on those nodes. *Calibration points used in determining divergence times, 
which corresponds to three nodes of phylogenetic tree in Fig. 1. 


tree were highly supported by ML bootstrapping and Baye- 
sian posterior probabilities (BPP). 

The ML tree showed that the two representative orders 
(Alcyonacea and Gorgonacea) of the subclass Octocorallia 
are not monophyletic group. The temperate soft coral D. 
gigantea was clustered with Pseudopterogorgia bipinnata 
(Gorgonacea: Gorgoniidae), instead of being in sister rela- 
tionship with Sarcophyton glaucum (Alcyonacea: Alcyonii- 
dae) (Fig. 1). High ML bootstrapping values and Bayesian 
posterior probabilities strongly suggest that molecular data 
do not support the conventional taxonomic system in this 
subclass that is mostly based on morphological similarity. 
Previous studies showed a similar incongruence between 
taxonomic system and molecular phylogeny particularly at 
an ordinal level and a lack of resolution at the base of octo- 
corallian phylogenetic tree. Recently McFadden et al. (2006) 
found that polyphyletic relationships are widely spread in 
the subclass Octocorallia, especially in five orders (Alcy- 
onacea, Gorgonacea, Stolonifera, Telestacea, and Protoal- 
cyonaria). Previously these orders were proposed to be 
combined in one order, Alcyonacea (Bayer, 1981) and this 
revision has been widely accepted (Fabricius, 2001). On the 
other hand a lack of resolution at deep nodes was interpreted 
as a possible result of radiating evolution among major line- 
ages of octocorallians (McFadden et al., 2006). The radiation 
hypothesis waits further examination and supporting by dif- 
ferent markers. Like the tree Medina et al. (2006) found, our 
tree reconfirmed that the corallimorpharians were nested 
within scleractinians. Based on this finding, it was argued 
that the taxonomy of scleractinian based on skeleton does 
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not define a natural group and scleractinian skeletons might 
be lost in the nonskeletonized corallimorpharian lineage 
(Medina et al., 2006). However this hypothesis was not 
upheld in other study by Fukami et al. (2008). Based on 
phylogenetic analyses of extensive taxon sampling and dif- 
ferent molecular data of nuclear genes (ribosomal DNA and 
B-tubulin) and mitochondrial genes (CYTB and COD), they 
found that the corallimorpharians and scleractinians are 
monophyletic group, respectively (Fukami et al., 2008). So 
this conflict warrants further investigations. For example, 
another phylogentic analysis using mtDNA genome data 
from the same sampling scheme used in Fukami et al.’s 
(2008) analysis might resolve the conflict, although not all 
the taxa are needed. Mitochondrial genome data from antho- 
zoans does not yet reach that extent. 

The ML tree topology in Fig. 1 is the same as Medina et 
al.’s (2006) except for the zoanthiniarian, Palythoa sp. 
According to Medina et al. (2006) this species was clustered 
to actiniarians with a low BPP value (84% and less then 50 
% in ML bootstrapping). Our result showed a relatively 
close relationship of zoanthiniarian to octocorallians than to 
actinarians, although the node still has low ML bootstrap- 
ping value (79%). It is likely that a difference in aligned 
data set for analyses between studies might be responsible 
for the different results. In addition, only two DNA genes 
are usable for Palythoa sp. (Table 1). So uncertainty on the 
position of Palythoa sp. still remains. 

We estimated divergence times of each node using the ML 
tree and the penalized likelihood method implemented in 
R8S program. As Medina et al. (2006) used two calibration 
schemes which were based on fossil records, we employed 
these two calibrations separately in the estimation for all the 
nodes except for the root and the next two nodes connected 
to the branch having low (79%) ML bootstrapping value. 
Ranges of divergence time for each node represent each re- 
sult from the two calibration schemes (Fig. 1, Table 3). The 
basal Astrangia node (70 Ma) in the “Robust/Short” clade 
consistently resulted in more elevated values than when the 
other nodes in “Complex/Long” clade were used. In overall 
as nodes go back to deep along the tree, the range becomes 
wider. Therefore a careful interpretation is needed for such 
values from deep nodes. Our estimates are slightly different 
from those of Medina et al. (2006). The difference mainly 
occurred around upper boundary of ranges (Table 3). Our 
rooting method is based on midpoint rooting which defines 
a root in the middle of the longest branch connecting any 
two terminal taxa. This was inevitable in our analysis be- 
cause of absence of relevant outgoup to the ingroup. On the 
contrary, Medina et al. (2006; see Fig. 1) used the octocoral- 
lians as outgroups to estimate diverging ages of two deeper 
nodes among the hexacorallians. Therefore it appears that 
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the midpoint rooting is related to our elevated estimates of 
divergence times. Molecular dating of divergence time has 
been known to have various sources of errors including 
phylogenetic uncertainty, substitutional noise, and lineage 
effects (Sanderson and Doyle, 2001). 

According to our analysis the octocorallians seem to start 
diverging around end of or after Cretaceous period (50-79 
Ma). This age relatively very early compared to the long 
history of stony sea corals (>240 Ma). However because 
the estimate is based on only four taxa of the Octocorallia, 
whether the octocorallians really diversified and radiated 
widely in such a recent past remains to be tested. Thus, it 
will be interesting how the result will be when octocoralian 
taxa are comprehensively examined. Since mitochondrial 
genome data from the subclass Octocorallia is very scarce 
up to date, the sequence information of D. gigantea can be 
used as an useful molecular and comparative framework for 
another mitochondrial genome sequencing and constructing 
a more comprehensive phylogeny of the soft coral, octoco- 
rallians. 
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